Factorization of Dephasing Process in Quantum Open System 
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The fluctuation-dissipation relation is well known for the quantum open system with energy 
dissipation. In this paper a similar underlying relation is found between the bath fluctuation and 
the dephasing of the quantum open system, of which energy is conserved, but the information is 
leaking into the bath. To obtain this relation we revisit the universal, but simple dephasing model 
with quantum non-demolition interaction between the bath and the open system. Then we show 
that the decoherence factor describing the dephasing process is factorized into two parts, to indicate 
the two sources of dephasing, the vacuum quantum fluctuation and the thermal excitations defined 
in the initial state of finite temperature. 

PACS numbers: 05.40.-a, 03.65.Yz, 32.80.-t, 42.50.-p 



I. INTRODUCTION 

A realistic quantum system can rarely be isolated from 
its surrounding environment (or called the "bath" ) com- 
pletely P, 0]. When it is coupled to the bath with a 
large number of degrees of freedom, decoherence hap- 
pens. There usually are two distinct decoherence effects 
of the bath on the quantum system: when the energy 
exchange is allowed by the interaction between the open 
system and the bath, the system energy usually dissi- 
pates into the environment irreversibly and we name this 
effect by quantum dissipation 0, 0, la] ; another effect is 
entitled by quantum dephasing, which occurs with no en- 
ergy exchange, but an irreversible process of information 
loss happens in the considered open system @, 0, II] ■ 

As to quantum dissipation, the well-known fluctuation- 
dissipation relation [J 0, E3, EH reveals to what extent 
the bath fluctuation depends on the quantum dissipation 
process of the open system, or how the phenomenological 
damping rate is determined microscopically by the ran- 
dom couplings to the degrees of freedom of the bath. As 
to quantum dephasing, however, there does not exist a 
similar statement clarified for the quantum decoherence 
only with dephasing and without loss of energy. 

This paper is devoted to find the intrinsic relation be- 
tween the pure dephasing and the some random nature of 
the bath. We begin with a general dephasing model for 
an open system interacting with a bath of many harmonic 
oscillators via a coupling of quantum non-demolition |T3 ] . 
This model is known to be universal in weak coupling 
limit Here, we characterize the dynamic process of 
dephsing with the so called decoherence factor, which 
linearly accompanies the off-diagonal elements of time 
evolution of the reduced density matrix for the system 
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obtained by tracing over the variables of the bath. We 
find that the decoherence factor is a product of two fac- 
tors, one of which is determined by the excitations of 
the bath while another does not vanish for the bath ini- 
tially in vacuum state. This observation clearly indicates 
that there exist two sources of quantum dephasing, which 
originate from both the vacuum quantum fluctuation and 
the thermal excitations of the heat bath respectively. 

The paper is organized as follows. In section II, we 
describe a universal model for quantum open system 
interacting with a bath of many bosons through non- 
demolition couplings. In section III, we show that a fac- 
torization structure appears in the dynamic dephasing 
process of the quantum open system. In section IV, we 
present the central result of this paper, i.e., the relation 
between dephasing and the fluctuation of the bath. In the 
section V, we study the dephasing of the quantum open 
system while the bath is in thermal equilibrium state. 
Finally, in the last section we conclude this paper with 
the remarks about the relationship between the quantum 
dephasing and the generalized thermalization due to the 
entanglements of the system with the bath. 



II. BOSON BATH FOR DEPHASING 

In general, quantum dephasing process is microscopi- 
cally considered as the vanishing of off-diagonal elements 
of the time evolution of reduced density matrix of an open 
system interacting with its surrounding environment or 
bath of infinite degrees of freedom. The simplest model is 
the composite system consisting of a system interacting 
with the bath of infinite bosons [H • 

The model Hamiltonian 

H = Hs + Hj + Hb 

is decomposed into three parts, the system part Hs, the 
bath part Hb = ^jO^aj with the bosonic creation 
(annihilation) operators a}- (ctj ) for the mode of frequen- 
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FIG. 1: The energy level configuration for the quantum 
open system defined within the subspace of excited states \n) 
spanned by the transition operators 0„ = |<7)(n| from the 
ground state \g). 



cies ujj (j — 1,2,...) and the interaction between the 
quantum open system and the bath [l3| 



H T = HG^i^aj +B..C.) 



(1) 



where we assume that h£j has the dimension of energy. 
The operator-volumed coupling G depends on the sys- 
tem variables. To conserve the energy of the open sys- 
tem, if dephasing can happen, the coupling is required 
to be of quantum non-demolition, i.e., [Hs,Hi] = and 
[Hb,Hi] ^ 0, or [Hs,G] = 0. In general the considered 
system has N energy levels with eigen-states \n) and cor- 
responding eigen-values Q n (n = l,2....N) [lj|, i.e., 

N 



H s = J2 m n\n)(r 

n=l 
N 

G = ^2g n \n)(n\, 



(2) 



n=l 



where g n is dimensionless since UG^j has the same di- 
mension of energy as that of K£j. The mode num- 
ber N may be infinite for a macroscopic heat bath. 
We assume that the system can be embedded into a 
N + 1— dimensional space with an extended basis vector 
\g) such that \n) = (ftAg) (see Fig. [1]) where the tran- 
sition operators <j) n — \g){n\ has commutation relations 
with the projection operators p n = \n)(n\ 



[<l>n,Pm] = S mn (j) n . 



(3) 



A typical example of this model is that the open system 
is a single mode boson with Hamiltonian Hs — ftu>owb 
defined by the bosonic creation (annihilation) operator 
W (b) for the mode of frequency uoq. The interaction 
Hi = htfbJ2j(£j a j+H- c -) between the boson and the 
bath [15( was frequently used to give a quantum approach 
for measurement of boson state. It originates from the 
generic oscillator coupling of the system coordinate q to 
the bath variables x,-, i.e., 



E to x i ~ E /; "- + ti bta l + ( 4 ) 

3 3 

in a large detuning limit \ujj — u)q\ 3> £j. 

We can exactly solve the Heisenberg equation for the 
operators <f> n and <Zj 

(j) n = -iQ n cj) n - iXg n cj) n , 

cij — —iujjCLj — iCjG (5) 

with the quantum noise operator 

X = Efe«i+ H.c). 

3 

Here, the noise operator satisfies the Brownian conditions 
in (X(t)) = and (X(t')X(t)) in an equilibrium 
state p and the thermodynamic average (...) is defined as 
{A) = Tr(pA). 

Since G is conservable or G (t) =G (0) = Go, we explic- 
itly obtain 



a j (t)=e- i ^ t a j (0)-iQr ]j (t)G 
where the phase operator in b(t) reads 



(6) 



^w = E^W a ^°)+ H - c -] 



and rjj (t) = i (e luJjt — l) fwj. Then the noise opera- 
tor can be expressed as a linear combination of initial 
operators Oj (0), (0) and Go, i.e., 



X(t)=Y(t)-G F(t) (7) 
where the "o-number" term 

nt) = Efc e "^(°)+ H - c -] ( g ) 



is due to the free evolution of the reservoir modes, and 
the c-number term 



. dw , „ / smut 
Fit) =2 / — J (to) (t 

UJ \ UJ 



(9) 



can be regarded as a time dependent "force" and defined 
by the spectral density function of the bath 



J(u;) = El0l 2 <^-^)' 



(10) 



This arises from the back action of the system on the 
bath. 
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III. FACTORIZING DEPHASING PROCESS 

To demonstrate the dynamic process of quantum de- 
phasing of the open system, we now calculate reduced 
density matrix for the time evolution of the open system. 
A pure dephasing process means that the off-diagonal el- 
ements of the reduced density matrix of the open system 
vanish, while the diagonal elements remain unchanged 
in such an ideal case. Usually the off-diagonal elements 
depend on both the initial state and the dynamic vari- 
ables of the bath and it can vanish in the thermodynamic 
limits. 

We assume that the initial state of the composite sys- 
tem is of the factorization form 

|V'(0)) = ^c n |n)(g)|{rn i }). 

n 

It means that the initial state of the open system is a 
superposition of the eigen states \n) = <jr n \g) while the 
bath is initially in Fock state ({m^}) = Yij \ m j)- Here 
rrij stands for the excitation number of Fock state \m,j). 
Let |0) = \g) £§) |{0j}), which is invariant under the op- 
eration of the evolution operator U (t) — exp(—iHt/?i), 
i.e., [/(t)|0> = |0>. 

Then, according to the explicit expressions ^ for op- 
erators <j) n and Oj, the time evolution of the composite 
system is calculated with a similar method in Ref. [16j | as, 

hK<)> = u(t)\i>(o)) 

= J2 Cn [U (t) ^t/t (t) ] u {t) \ g) g, \ {m , }) 
n 

= ^c^t (-t)| ff >®|{m i }>e- i ^W 

71 

= ^^(t) |n> »e i5 " Z( -' ) IK}) 

n 

or 

\i>(t))=J2 c n(t)\n)®\l3 n ). (11) 

n 

Here, 

c„ (t) = Cne -^t e iglF(t) e -iY lim ^ j t 
and the coherent state 

|/3n}=n^^na i )K) (12) 

3 

is defined by the displacement operators Dj (a) — 
exp(aaj — a*a,j) with 

«i = (t) = ^ (e-^ 4 - 1) ■ 

From the time evolution of density matrix p (t) = 
\ip (t)) (ip (t)\ for the composite system, we calculate the 



reduced density matrix of the open system 

Pa (*) = X! C « C ™ l n ^"l 

71 

+ J2 CnCe i9 -W (P m \p n ) \n) (m| (13) 

by tracing over the variables of the bath. Here 

mn (t) = (fl m -(l n )t+ {g 2 n - g 2 m ) F (f) 

is the time dependent real number. 

To characterize the coherence of the quantum open 
system, we use the decoherence factor [f| 

D n , m (t) = (/3 m |/3„) = Yl D mn [ mj ]. (14) 
3 

Initially, the decoherence factor D„ jTn (0) = 1, i.e., the 
system owns a completely ideal quantum coherence. 
Here, each of two factors in the decoherence factor 

D m n[m ] = {m \D ((g n - g m ) ay) \m d ) , 

is calculated as 

D mn [TOj] = (T 3 Z ™"^ (4) L m] (z nm .j(t) ) (15) 

in terms of the Laguerre polynomial L m , (z nm -j(t)) of 
variable 

Znm;j(t) = (flVi - 9™f \$jVj 

With the above results, we can observe that the deco- 
herence factor can be decomposed into two parts, i.e., 

D n , m (t) = D^l(t)D^(t) (16) 

where 

^m(*)=n e " lwt) 

3 

originates from the vacuum quantum fluctuation of the 
bath and 

D^(t)=l[L mj (z nm . tj (t)) (17) 

3 

means the bath excitation. Actually, when the bath 
is initially in the vacuum state \{nrij = 0}), the sec- 
ond factor D^m 3 ^ (t) becomes a unity operator since 
L mj=0 (x) = 1. The decoherence factor D n ^ m (t) = 

Dn)m (t) means that the dephasing only results from the 
vacuum fluctuation of the bath without thermal excita- 
tion. When the initial state of the bath (\{rrij}}) is oc- 
cupied by large amount of excitation, the macroscopic 
feature of the bath (high excitations) can induce the de- 
phasing of the open system. This is because L m . (x) 
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approaches the zero-order Bessel function Jo (%) when 
irij — » oo 17]. Then 

Dmn K] = e -5 2 — «W J (Wiif*)) ■ ( 18 ) 

Since the Bessel function with real variables is a decay- 
ing oscillating function, the decoherence factor D n m (t) 
(|16p approaches zero when t tends to infinity. We will 
discuss how a thermal equilibrium state can induce the 
fast dephasing. The above arguments show that, even 
though the system energy is conserved, there still exist 
the two sources of the pure quantum dephasing induced 
by the bath, i.e., the vacuum quantum fluctuation and 
the thermal excitations, which leak the information of 
the system into the bath. 

We can explicitly demonstrate these results in the con- 
tinuous limit. As a illustration we take Ohmic spectral 
density of the bath [H 

J (w) = 7we-" /r (19) 

to make the sum in 

D% (*) = e-( ff »-^) 2 W%MI 2 /2 (20) 
as a integral 

El^WI 2 = / ^H^^y- ( 21 ) 

3 

Here, 7 is a dimensionless coupling constant, and T the 
bath's response frequency. When the initial state of the 
bath is in vacuum state, the decoherence factor in Eq. (fl6|) 
becomes 

A*, m (*) = (l + r¥^ ( " 3 \ (22) 

Figure [2] shows that the decoherence factor (|22|) , the de- 
cay of which is similar to an exponential decay. Actually, 
at short time limit, T 2 t 2 <C 1, the decoherence factor (|22|) 
becomes Gaussian decaying faction as 

D n , m (t) = e -W. (23) 

Here the characteristic time r of the dephasing of the 
open system is defined by 

t~ x = On -fa)rj-. 



IV. DEPHASING-FLUCTUATION RELATION 

Now we take another approach to calculate the deco- 
herence factor, which will offer us a new angle to un- 
derstand the source of dephasing. Here, we adopt the 

notation of the standard deviation A A = J (A 2 ) — (A) 2 




0. 5 1 1.5 2 

t 



FIG. 2: Plot of D„, m (i) (decoherence factor) (l22l) vs t (time) 
for various value of the coupling parameter 7. \g m — g n \ = 1 
and r = 1. 

for a given operator A. We write down the decoherence 
factor 

D n , m (t) = ({I j }\e i +V\{I j }) (24) 

in terms of the phase difference operator 

<M*) = (5™ -9m)Z(-t) 

for the bath initially with a factorized state \{Ij}) = 
Hj We denote the average ({/,•} \A\ {/,}) by (A). 

For the small variation \<f> (t) — (</>(£)} | of the phase <j> (t) 
around it average {4>(t)) over any state we have approxi- 
mately [l|] 

D n , m (t) = (e^W) ~ e ^*(*)> e -i( A ^)) 2 . (25) 

It shows that the phase fluctuation A(j>(t) results in the 
loss of quantum coherence or quantum dephasing char- 
acterized by the decaying of the decoherence factor. For 
the initial Fock state with |J,-) = \m,j) we discussed in 
the last section, (</>(t)) = and the phase fluctuation due 
to the bath fluctuation can be separated explicitly into 
two parts, i.e., 

(A0(t)) 2 = (A4>(t))l + (A<Hi))/ ■ (26) 
where the vacuum fluctuation part is 

(A0W)o = E z «^W ( 2? ) 

3 

and the bath excitation part 

(A0(i))J = J2 2m 3 Z nm;j(t). (28) 
3 

Actually the result is the same as what we have ob- 
tained above in some cases. This is because the cou- 
plings are linear with respect to the operators Oj and a] . 
Through some simple calculations, we obtain the deco- 
herence factor as 

D n . m (t) = e -h(^(t))l e -¥^mf f 

= D^l{t)D^ m {t) (29) 
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where the thermal excitation part can be rewritten as 



(30) 



When the bath is initially in vacuum state \{rrij}) — 
\{rrij =0}), fj (rrij) — 0. Then the decoherence fac- 
tor becomes D n;m (t) = Dn ln (t). Otherwise we have 
fj (rrij) > 0, it means that the decoherence factor in 
Eq. pT))) will decay with time . 

Next we consider the difference of the decoherence fac- 
tors obtained through different approaches in Eq. (fTB|) and 
Eq. ([29|) respectively. When rrij and (<?„ — g m ) /ujj are 
all small, i.e., the bath is setup in a low-excitation state 
and coupling is weak, the Laguerre polynomial in Eq. (117p 
can be approximately as 

L rnj (x) a e~ m ^. (31) 

Then the two decoherence factor in Eq. p^)) and Eq. p5|) 
has the same form. We consider the time evolution of 
the average excitation number of the bath, i.e., 

N b (t) = J2 (°i (*) a i (*)) = Nb (°) + 5Nb (*) ( 32 ) 

3 

where the initial average excitation number of the bath 
is N B (0) = J2j m j> and 

3 

= (G 2 ) Jdu;J(u,)±sm*f (33) 

is the quantum fluctuation of the bath excitation num- 
bers, which is independent of the initial state of the bath. 
In the above equation, the term (Gq) indicates the aver- 
age of the square of the population number of the open 
system. This means 



£&»&(*)i a = 



SN B (t) 
(G 2 ) ■ 



(34) 



Then the decoherence factor (1161) can be rewritten in 
terms of the fluctuation of excitation number, 



D n , m (t) oc exp 



1 



(<7n ~ 9m) 



2 SN B (t) 
(G 2 ) 



(35) 



This equation shows that the decoherence factor is de- 
termined by the quantum fluctuation of bath 8N B (t) and 
(Gq). When the initial state of the open system is given, 
the value of (Gq) is fixed in the mean while. That is to 
say, (Gq) is not relevant to the source of dephasing of 
the open system, but the energy of the open system will 
contribute to the dephasing. Thus we conclude that the 
quantum fluctuation of the bath excitation characterized 
by SN B (t) induces the dephasing of the open system. 



V. DEPHASING IN THERMAL EQUILIBRIUM 
STATE 



When the bath is prepared in a pure state we have 
found that the dephasing-fluctuation relation for quan- 
tum dephasing is similar to the well-known fluctuation- 
dissipation relation for quantum dissipation. Corre- 
spondingly the dephasing process can be understood sep- 
arately according to the quantum fluctuation and the 
pure state excitation of the bath. In this section we will 
show that these observations also hold exactly for the 
case that the bath is initially in a thermal equilibrium 
state. 

For the dephasing problem of the open system at finite 
temperature, the thermal equilibrium bath is described 
by the density matrix 



Pb = 



n 



(36) 



Initially, the density operator of the composite system is 
a direct product 

p(0) = hK0))(V (0) | 8) p B 

with the initial state \ip (0)) = ^2 n c n \n) of the open sys- 
tem. Using the coherent-state representation, the density 
operator can be rewritten as 



with the P— representation for diagonal elements 

1 



P(A;) = 



Tt(mj) 



a -|A 3 | /(m 3 > 



(37) 



(38) 



-l . 



where (rrij) = ^ nrriiUJj — l) "is the average excitation 
number in the mode of the frequency LOj. 

By applying the results obtained in the previous sec- 
tions, we can calculate the reduced density matrix 

Ps (t) = C » C « I") H+ E C " W °*m (*) I") H D nl (t) ■ 

n n^m 

(39) 

to characterize the quantum coherence of the open sys- 
tem. Then the decoherence factor in the above equation 
is factorized as [13 



4Tiw=n/ rf2 ^ 



TT{rrij) 

J 

=^(*)n e " wcn 

3 

where a jmn = (g n - g m ) a 3 and 
9j (T) 



(\j\ Dj (a jmn ) \Xj) 
(40) 

(41) 
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FIG. 3: Plot of D l n,L (t) (decoherence factor) igOj vs t (time) 
and T (temperature). \g m — g n \ = 1,7 = 1 and F = F 

Figure[3]shows that the decoherence factor D n . m (t) de- 
pends on the time t and the temperature T of the bath. 
When the temperature T increases, the decaying deco- 
herence factor will be enhanced, i.e., the open system 
loses its quantum coherence. 

When the temperature approaches the absolute 
zero degree {[3 — > oo), we approximately have 
^phujj _ 1 ) 1 ~ e -/3?iwj Then the decoherence factor 
PO)) becomes 

D^ m (t) = DS? M (t) exp (- ^-^/^ klT*] . 

(42) 

It shows that, in the low-temperature limit, the above 
decoherence factor exponentially decays as the square of 
temperature T 2 increases, which is called Gaussian de- 
cay. When the bath is prepared in the high temperature 
(Jihujj — > 0), we have e° i — 1 ~ fihujj, and then 

5i (T) = ( P „- 5m ) 2 %^ocT. (43) 

It means that, when the bath approaches the high- 
temperature limit, the decoherence factor (|40[) exponen- 
tially decays as the temperature T increases. 




0.2- 



0.0-| , , --; ■, r j ■, ; ^"""r^v"""l 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 
t 

FIG. 4: Plot of D [ J} m (t) (decoherence factor) igDJ vs t (time) 
for various value of temperature T. \g m — g n \ = 1, 7 = 1 and 
F = 1. 



In addition, Figure[4]shows that the rise of temperature 
will accelerate dephasing described by Ea. (j4H)) . 



VI. CONCLUSION WITH REMARKS 

In summary, through a universal model for a quan- 
tum open system with non-demolition coupling to the 
bath, we find the intrinsic relation between the dephas- 
ing of the open system and the quantum and thermal 
fluctuation of its bath. Usually, the couplings of a quan- 
tum open system to its environment is very complicated, 
and then intuitively the dephasing process should depend 
on the details of system-bath couplings. So generally 
the present model used in this paper seems to be too 
oversimplified. However, we would like to mention the 
Caldeira and Leggett's studies about quantum dissipa- 
tion [3j , which proves that any bath with weak couplings 
to the open system can be modeled as a collection of 
non-interacting bosons. As for the quantum dephasing, 
we even found a similar conclusion to universally consider 
the quantum dephasing problem in quantum computing 
[2fj| . Thus the present investigation can also be consid- 
ered to be universal, and dephasing process is generally 
clarified as the two kinds of origins, the quantum fluctu- 
ation in the vacuum and the thermal excitation in finite 
temperature. 

Before concluding this paper, it is necessary to say 
some words about the repletion and difference between 
the dynamic dephasing process studied in this paper and 
the thermalization to the thermal equilibrium induced 
by the bath [U [H, SI H| . First, we have to note 
that the dephasing process in this paper can describe the 
thermal equilibrium with canonical state in a straight- 
forward way. Most recently people revisit the investi- 
gation to explore the possibility replacing the equal a 
prior probability postulate in statistical mechanics by a 
general canonical principle (2lj |: an arbitrary entangling 
pure state of the total system consisting of the "small" 
system plus "large bath" can be traced over the variables 
of the bath to give a generalized canonical state due to 
the "Large Number Law" or overwhelming majority of 
wave functions in the subspace by obeying some global 
constraint for the "universe" . Actually, if this global con- 
straint is particularly determined by the energy interval 
encompassed by the microcanonical ensemble, the above 
mentioned "tracing" operation just results in the ther- 
mal equilibrium distribution. Such thermalization pro- 
cess can also be described as a typical quantum dissipa- 
tion phenomenon, a dynamically quantum process with 
energy exchange between the system and the bath 
However, there is not energy exchange in our present de- 
phasing model and thus it give rise to a mixture rather 
than the thermal equilibrium distribution. Maybe there 
exist another type global constraint other than the en- 
ergy interval, we believe that it is an open question, which 
could be solved in a general framework. We will continue 
the exploration along this line in future studies. 
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